Spiral ground state in the quasi-two-dimensional spin-| system Cu2Ge04 



O 

(N 



Alexander A. Tsirlin/'Q Ronald Zinke, 2 Johannes Richter, 2 and Helge Rosner 1 '!]] 

1 Max Planck Institute for Chemical Physics of Solids, Nothnitzer Str. 40, 01187 Dresden, Germany 
2 Institute for Theoretical Physics, University of Magdeburg, P.O. Box 4120, 39016 Magdeburg, Germany 

We apply density functional theory band structure calculations, the coupled-cluster method, 
and exact diagonalization to investigate the microscopic magnetic model of the spin-i compound 
Cu2Ge04. The model is quasi- two-dimensional, with uniform spin chains along one direction and 
frustrated spin chains along the other direction. The coupling along the uniform chains is antifer- 
romagnetic, J ~ 130 K. The couplings along the frustrated chains are Ji ~ —60 K and J2 — 80 K 
between nearest neighbors and next-nearest neighbors, respectively. The ground state of the quan- 
tum model is a spiral, with the reduced sublattice magnetization of 0.62 fiB and the pitch angle of 
84°, both renormalized by quantum effects. The proposed spiral ground state of Cu2Ge04 opens a 
way to magnetoelectric effects in this compound. 

PACS numbers: 75.30.Et, 75.10.Jm, 71.20.Ps, 75.50.Ee 
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I. INTRODUCTION 

Quantum magnetism is a field of fundamental research 
focused on exotic ground states and non-trivial low- 
temperature properties ^ Nevertheless, certain effects 
in quantum magnets are also relevant for applications. 
Spin-chain compounds show ballistic regime of heat 
transport 4 whereas frustrated magnets are capable of a 
strong magnetocaloric effect 4^ Additionally, many of the 
frustrated magnets undergo spiral or, in general, incom- 
mensurate ordering, and reveal ferroelectricity induced 
by a magnetic fields A frustrated spin chain with com- 
peting ferromagnetic (FM) nearest-neighbor (Ji) and an- 
tiferromagnetic (AFM) next-nearest-neighbor (J 2) cou- 
plings is the simplest spin model giving rise to spiral mag- 
netic correlations at J2/J1 < — \ (Ref. 0). This model is 
easily realized experimentally and has a clear structural 
footprint, a chain of edge-sharing CuX 4 plaquettes with 
X being oxygen^r— chlorine^ 2 - or even nitrogenJ^ Such 
chains typically show FM J\ due to the nearly 90° Cu-X- 
Cu angle and AFM J2 due to the Cu-X-X-Cu superex- 
change. Indeed, many compounds of this type undergo 
spiral magnetic ordering and sometimes exhibit magnetic 
field-induced ferroelectricityJ4 However, the detailed mi- 
croscopic understanding of these effects remains challeng- 
ing, and even the electronic origin of ferroelectricity in 
spin-chain cuprates is vividly debated^ 

Interchain couplings are an important feature of any 
real material. The couplings between spin chains can 
modify the ground state qualitatively by inducing a long- 
range order with finite sublattice magnetization! 16 i 17 In 
the case of frustrated spin chains, such couplings influ- 
ence the behavior of doped systems^ and play a decisive 
role for the stability of exotic phases in high magnetic 
fields i 19 ' 20 Regarding magnetoelectric effects, the inter- 
chain couplings naturally determine their temperature 
scale by adjusting the magnetic ordering temperature. 

Theoretical studies of coupled frustrated spin chains 
remain a challenge owing to the two-dimensional (2D) 
and frustrated nature of the problem. Therefore, experi- 



mental benchmarks are especially important. The avail- 
able frustrated-spin-chain compounds show relatively 
weak interchain couplings,— while the relevance of the 
opposite regime with strongly coupled frustrated spin 
chains remains unclear. A common and a somewhat 
naive picture suggests that leading exchange couplings 
should run along the structural chains owing to shorter 
Cu-Cu distances.— 

In the following, we present a microscopic magnetic 
model of Cu 2 Gc04. This compound is a unique example 
of a 2D system of strongly coupled frustrated spin chains. 
The coupling J between the frustrated chains is so strong 
that the system can be equally viewed as uniform spin 
chains along J with the frustrated interchain couplings J\ 
and J 2 (see Fig. [T]). Both descriptions relate to certain 
features of the magnetic behavior: while the uniform- 
chain model fits the magnetic susceptibility of CU2GCO4 
down to T/J ~ 0.5, the ground state of the 2D model is 
a spiral, which is typical for the frustrated Ji — J2 spin 
chains. 

The crystal structure of Cu2Ge0 4 belongs to the 
spinel typei^ 2 - Magnetic properties were studied in a rela- 
tion to the spin-Peierls compound CuGe03. The low- 
dimensional magnetic behavior of Cu2Ge0 4 resembles 
CuGe03 indeed. However, no signatures of the struc- 
tural distortion or spin gap were found down to 10 K, 
and the long-range magnetic ordering at Tjv = 33.1 K is 
observed instead^ 

Yamada et al£2- analyzed Cu 2 Ge04 using the 
anisotropic pyrochlore lattice model with two inequiv- 
alent exchange couplings that are J\ and J c in our nota- 
tion (upper left panel of Fig. [IJ. This model arises from 
a straight-forward and naive geometrical consideration of 
the spinel structure, with inequivalent couplings driven 
by the tetragonal distortion of the parent cubic system. 
At J c /J\ -C 1, the anisotropic py rochlore lattice splits 
into chains. According to Ref. |23, Cu2Ge0 4 is close to 
this limit, with J\ = 135 K and J c / J\ = 0.16. Starykh 
et al^ studied the 2D analog of the model theoretically, 
and proposed a quantum-disordered valence-bond-solid 
ground state. 
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FIG. 1. (Color online) Top panel: crystal structure of 
Cu2Ge04 (left) and a single magnetic layer in the ab plane 
(right). Bottom panel: a sketch of the spin spiral with the 
pitch angle 7 (left), and the magnetic model of Ji — J2 frus- 
trated spin chains coupled by J (right). Circles and dots 
denote the positions of the Cu atoms. Lines in the top left 
panel show the anisotropic pyrochlore lattice considered in 
Ref . HI 



II. BAND STRUCTURE 

As a derivative of the spinel structure, Cu2GeC>4 might 
be thought of as a three-dimensional network of CuOg 
octahedra. However, this description ignores essen- 
tial features of the electronic structure. In oxide com- 
pounds, Cu+ 2 tends to adopt a four-fold coordination 
(CuC>4 plaquette) having dramatic influence on the or- 
bital ground state and magnetic properties. Such pla- 
quettes can be recognized in Cu2Ge04, and lead to a 
peculiar supcrexchange scenario. Four short bonds to 
oxygen (1.95 A) form the CUO4 plaqucttcs in the ab 
plane, whereas the two remaining Cu-0 bonds are much 
longer (2.50 A). Edge-sharing CUO4 plaquettes comprise 
structural chains that run along a or 6, with parallel 
chains forming layers in the ab plane (upper right panel 
of Fig. [T|). Equivalent layers with differently directed 
structural chains alternate along the c axis. In the fol- 
lowing, we denote the direction of the structural chains 
as b' and the perpendicular direction as a', to distinguish 
those from the crystallographic a and b axes. GCO4 tetra- 
hedra connect the chains into a three-dimensional (3D) 
framework (Fig. [1]). 

To evaluate individual exchange couplings, we perform 
scalar-relativistic density functional theory (DFT) band 
structure calculations using the FPL09 . 00-33 code^ We 
apply the local density approximation (LDA) with the 
exchange-correlation potential by Perdew and Wang^ 
and use a well-converged k mesh comprising 3350 points 
in the symmetry-irreducible part of the first Brillouin 
zone. With LDA calculations, we are able to identify 
relevant states, and to evaluate hopping parameters ti 
via a fit with an effective one-orbital tight-binding (TB) 
model. The hopping parameters are introduced into a 



TABLE I. Leading exchange couplings in Cu2Ge04: hopping 
parameters ti of the TB model, AFM contributions to the 
exchange couplings J/^ FM = 4t;/[/eff, and the total exchange 
integrals Ji from LSDA+f/ calculations with Ud = 6.5 eV. 
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Hubbard model with the effective on-site Coulomb repul- 
sion potential U c g = 4.5 eV.— ~— In the case of low-lying 
excitations, the Hubbard model is further reduced to a 
Heisenberg model under the conditions of half-filling and 
strong correlations (ij -C U c g). Then, the AFM parts of 
the exchange integrals are evaluated as J^ FM — At^/U e g. 

An alternative way to evaluate the exchange couplings 
is to treat the strong correlations within DFT, via the 
mean-field-like LSDA+[/ approach. We calculate total 
energies for a set of collinear spin configurations, and 
map these energies onto a classical Heisenberg model. 
Thus, total exchange integrals Ji are estimated. In the 
LSDA+J7 calculations, we use the Coulomb repulsion 
and exchange parameters Ud = 6.5 ± 1 eV and Jd = 1 eV, 
respectively^ - — The double-counting-correction (DCC) 
scheme was set to the around-mean-field (AMF) option. 
The application of the fully-localized-limit (FLL) DCC 
had little effect on the exchange couplings. 

The LDA energy spectrum of Cu2Ge04 is typical 
for Cu +2 oxides. The mixed Cu 3d - O 2p valence 
bands extend down to —8 eV (Fig. [5]), with the states 
near the Fermi level predominantly formed by the Cu 
d x 2_ y 2 orbital (here, x and y align with the short Cu- 
O bonds). Germanium orbitals contribute to the bands 
around —10 eV, and show negligible DOS at higher ener- 
gies. While LDA yields a metallic energy spectrum due to 
the underestimation of electronic correlations in the Cu 
3d shell, LSDA+C/ restores the insulating scenario with 
the band gap of E g = 2.0 ± 0.3 eV for U d = 6.5 ± 1 eV. 

The Cu d x 2_ y 2 states are represented by four bands 
crossing the Fermi level and arising from four Cu atoms in 
the primitive cell of Cu2Ge04 (Fig. [3]). These bands are 
separated from the rest of the valence bands by a pseudo- 
gap. To extract hopping parameters, we construct Wan- 
nier functions (WFs) based on the Cu d x 2_ y 2 character M- 
This analysis evidences sizable nearest- neighbor (ti) and 
next-nearest-neighbor fa) hoppings along the structural 
chains. However, the hopping t along a' is comparable to 
ti and i2- Additionally, a weak diagonal hopping in the 
ab plane is found (Table HJ. The nearest-neighbor hop- 
pings perpendicular to the ab plane (t c ) arc —11 meV, 
yielding J^ FM as low as 1 K. The weak dispersion of 
the bands along T — Z also shows the pronounced two- 
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FIG. 2. (Color online) LDA density of states for Cu 2 Ge04. 
The Fermi level is at zero energy. 



dimensionality of the system. Introducing the hoppings 
into an effective one-band Hubbard model, we evaluate 
AFM parts of the exchange integrals J AFM (Table HJ). 

LSDA+C7 calculations modify the LDA-based scenario. 
We find FM nearest-neighbor coupling within the struc- 
tural chains, J x = -60 T 10 K for U d = 6.5 ± 1 eV. 
The next-nearest-neighbor intra-chain coupling J 2 = 
80 T 20 K and the interchain coupling J = 130 T- 30 K are 
basically unchanged. Further couplings in the ab plane 
are below 10 K. The interplane coupling becomes FM and 
remains weak. Thus, we establish the quasi-2D J— Ji— J 2 
model with a weak interlayer coupling J c (Fig.[l|. 

The quasi-2D model of Cu 2 Ge04 results from the 
strong tetragonal distortion of the spinel structure. The 
plaqucttc description (Fig. [I}, with the magnetic d x 2_ y 2 
orbital coplanar to the CUO4 plaquette, clarifies the 2D 
nature of the system. The couplings J c connect the pla- 
quettes lying in different planes, and therefore remain 
weak. By contrast, three sizable couplings in the ab plane 
establish a frustrated spin lattice. Our model is dissimilar 
to the anisotropic pyrochlore lattice proposed by Yamada 
et al££ The pyrochlore spin lattice omits the relevant ex- 
changes J and J 2 , and should be discarded. Cu 2 Ge04 is 
a frustrated magnet indeed, but the strong frustration is 
found in the J\ — J 2 chains rather than tetrahedral units. 

The FM nearest-neighbor coupling J\ should be re- 
ferred to the Cu-O-Cu angle of 91.8°. The microscopic 
origin of ferromagnetism is the Hund's coupling on the 
oxygen site£ The next-nearest-neighbor coupling J 2 is 
the AFM Cu-O-O-Cu superexchange. Similar values 
of 50 — 100 K for \J\\ and J 2 have been established for 
the archetype frustratcd-spin-chain compounds, such as 
LiCu 2 2 and LiCuV0 4 ^ 

Another remark on the structural implementation of 
the spin model regards the origin of the long-range cou- 
plings J and J 2 . Since the Ge orbitals weakly contribute 
to the valence states, both couplings should be assigned 
to a Cu-O-O-Cu superexchange. Despite an identical 
Cu-Cu distance (Table H]), a larger J value is caused by 
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FIG. 3. (Color online) LDA band structure of Cu2Ge04 (thin 
light lines) and the fit of the TB model (thick dark lines). The 
Fermi level is at zero energy. The fc-path is defined as fol- 
lows: T(0, 0, 0), X(0.5, 0, 0), A/(0.25, 0.25, 0), and Z(0, 0, 0.5), 
where the coordinates are given in units of the reciprocal lat- 
tice parameters 4ir /a and 4ir/c. 



the co-planar arrangement of the plaquettes in the adja- 
cent chains. By contrast, the next-nearest-neighbor pla- 
quettes within the chain ( J 2 ) lie in different planes due to 
the buckled chain geometry (Fig. [T]) . It is worth to note 
that the ab projections of the Cu 2 Ge04 and LiCuVC>4 
structures are very similar. However, LiCuV04 is a 
quasi- ID system with J <C |Ji|, J 2 , while the spin sys- 
tem of Cu 2 Ge04 is quasi-2D^ 



III. MICROSCOPIC MODEL 

In the following, we explore the ground state and finite- 
temperature properties of our model. We first consider 
the purely 2D regime described by the Hamiltonian: 

n i 

+ / ] J Sj,n ■ Sj, ra+ 1, (1) 
i n 

where the index n labels the structural chains (along &'), 
and i denotes the lattice sites within a chain n. The effect 
of the interlayer coupling J c is discussed in Scc. lIIIDl 

Our model can be viewed as frustrated J\ — J 2 chains 
(along b') which are uniformly coupled by J (along a'). 
Alternatively, one finds uniform spin chains along a' with 
frustrated interchain couplings J\ and J 2 along b' . While 
any of the parent ID models is rather easy to handle, 
a rigorous treatment of their 2D combination is a chal- 
lenging problem. Below, we apply the Lanczos diagonal- 
ization and coupled cluster method to achieve an accu- 
rate description of the ground state. By contrast, finite- 
temperature properties of the quantum model can only 
be accessed at high temperatures by a series expansion 
(HTSE), whereas conventional techniques, such as quan- 
tum Monte-Carlo or exact diagonalization, fail because 
of the sign problem or finite-size effects^ 
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A. Magnetic susceptibility 

Since the experimental information on Cu 2 Ge04 is 
restricted to the magnetic susceptibility and heat ca- 
pacity measurements in Ref. l23l . we discuss thermody- 
namic properties first. The experimental specific heat 
contains an unknown phonon contribution, hence the 
magnetic part can not be separated. Therefore, the 
magnetic susceptibility x(T) (Fig. [4]) remains the only 
quantity suitable for the comparison between theory and 
experiment. The estimated Curie- Weiss temperature 
9 ~ i(J + Ji + J 2 ) = 75 K is in good agreement with 
the experimental value of = 89 K 32 

For a further comparison, we derive the HTSE for our 
model: 

= N a9 2 ij% f J + J1+ J 2 J 2 + Jf + Jl \~ 1 
X 4k B T \ 2T AT 2 J ' 

(2) 

where Na is Avogadro's number, hb is Bohr magneton, a 
is the g-factor, and we used the expressions from Ref. |33| 
up to the third order in temperature. The calculated ex- 
change couplings (Table IJ) are in reasonable agreement 
with the experimental data down to 150 K (Fig. |4]). The 
deviations at lower temperatures are likely related to the 
divergence of the HTSE at T < J. To improve the fit 
at higher temperatures, a slight adjustment of the ex- 
change couplings is required. However, the third-order 
HTSE contains two J-dcpcndent terms only, hence an 
unconstrained fit of three exchange parameters is impos- 
sible. 

To access temperatures below 150 K, a simplification 
of the model is required. Since J exceeds \ J\\ and J2, the 
spin lattice is, to a first approximation, a set of uniform 
spin chains along a' . The respective ID model^ fits the 
experimental magnetic susceptibility down to 70 K with 
J = 140 K and g = 2.22. A similar fit with J = 135 K 
has been given in Ref. 0. However, Yamada et aL— er- 
roneously assign the spin chains to the structural chains 
(in their notation, J corresponds to J-y). Our DFT calcu- 
lations show that the uniform spin chains run perpendic- 
ular to the structural chains, whereas J\ is FM. Such an 
intricate situation is not uncommon for low-dimensional 
magnets, see Refs. EU andEU for similar examples. 

Below 70 K, the ID uniform-chain model overestimates 
the magnetic susceptibility of Cu2Ge04. This feature in- 
dicates an onset of 2D spin correlations. In contrast to 
the uniform spin-i chain having finite susceptibility at 
zero temperature, 2D and 3D systems usually develop 
a long-range order with vanishing susceptibility at low 
temperatures. The onset temperature of 2D spin correla- 
tions is a rough measure of interchain couplings. Indeed, 
the temperature of 70 K conforms to our estimates of 
Ji = —60 K and J2 = 80 K, the couplings between the 
uniform spin chains. 
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FIG. 4. (Color online) Fit of the experimental magnetic sus- 
ceptibility data with the uniform chain model (dashed line) 
and the comparison to the HTSE of Eq. @ (solid line) . Ex- 
perimental data are from Ref. [23|. 

B. Coupled cluster method 

The coupled cluster method (CCM) and its application 
to frustrated spin systems have be en p reviously reviewed 
in several articles, see, e.g., Refs. El IMS Therefore, 
we will give only a brief illustration of the main relevant 
features of the method. For more general information on 
the methodology of the CCM, see, e.g., Refs. |M andHU, 
and references therein. 

The CCM is a universal quantum many-body method. 
The starting point for a CCM calculation is the choice 
of a normalized reference or model state |<j>), together 
with a complete set of (mutually commuting) multi- 
configurational creation operators {C^ } and the corre- 
sponding set of their Hermitian adjoints {Cj,}. The CCM 
parametrizations of the ket- and bra- GSs are given by 

|*}=e s |$), S = J2SiCl; 
<*| = (<t>\Se- s , S=1 + ^2SiCi. (3) 

Using ($|C+ = = C/|$> VI ^ 0, C$ = 1, the com- 
mutation rules [Cl,C~k] = = [Cl,Ck], the orthonor- 
mality condition (&\CiCj = Sjj, and completeness 

Y^C+mmd = 1 = |*)(*| +£C+|*)(*|C7, we get 

a set of non- linear and linear equations for the correlation 
coefficients Si and Si, respectively. We choose a refer- 
ence state corresponding to the classical state of the spin 
model, i.e., a non-collinear reference state with up-down 
Neel -type correlations along the a'-dircction (uniform J 
chains) and with spiral correlations along the ^'-direction 
(frustrated J\ — J 2 chains). The spiral correlations are 
characterized by a pitch angle 7, i.e. |<E>) = |$(7)). In the 
quantum model, the pitch angle is typically different from 
the corresponding classical value 7 c i. Hence, we do not 
choose the classical result for the pitch angle, and rather 



5 



consider 7 as a free parameter in the CCM calculation. 
The value of 7 has to be determined by the minimization 
of the GS energy (in a certain CCM approximation, see 
below) given by £(7) = (3>(7)|e _,S -ffe' s, |<I>(7)), i.e., from 
the dE I ^7| 7=7qu = condition. 

In order to find an appropriate set of creation opera- 
tors, it is convenient to perform a rotation of the local 
axes on each of the spins so that all spins in the reference 
state align with the negative ^-direction. This rotation 
by an appropriate local angle 8i, n = #i in (7) of the spin 
on the lattice site (i,n) is equivalent to the spin-operator 
transformation 



cos 8i 



- sin <5i,„sf ; „ + cos 8i. n s 



(4) 



The reference state and the corresponding creation oper- 
ators are given by 



°L — S i,n 



s + s + 



|$> = UUI---> 

(5) 

where the indices (i, n), (j, to), (fc, I), . . . denote arbitrary 
lattice sites. This specified form of the creation opera- 
tors C£ and the corresponding reference state |$) imme- 
diately make clear that the general relations listed below 
Eq. @ are fulfilled. In the rotated coordinate frame, the 
Heisenberg Hamiltonian acquires a dependence on the 
pitch angle 7 (see Ref. [2l| for more details). 

The order parameter (sublattice magnetization) in 
the rotated coordinate frame is given by to = 

— l/iVXif n {*l s ? nl*)- The only approximation of the 
CCM is the truncation of the expansion of the corre- 
lation operators S and S. We use the well-established 
LSUBn scheme, where all multispin correlations on the 
lattice with n or fewer contiguous sites are taken into 
account. 

In contrast to Ref. [2l], which is focused on the J < 
Ji I , J2 regime, we consider the case of J > Ji | , J 2 . We 
also evaluate LSUBn approximations of higher order, up 
to n = 8. In the highest order of approximation, LSUB8, 
we have 21124 configurations, i.e., 21124 coupled non- 
linear equations have to be solved numerically. More- 
over, the minimum of £'(7) has to be found numerically 
to determine the quantum pitch angle 7 qu . For the nu- 
merical calculations, we use the program package CCCM 
by D.J.J. Farnell and J. Schulenburg42, 

Since the LSUBn becomes exact for n — > 00, the 
numerical result can be improved by extrapolating the 
"raw" LSUBn data to n — > 00 using the expression 
m n = moo + a/n + b/n 2 , cf. Refs. HI, Ell, EE andllsl 



C. Ground state 

To find the classical ground state of the model given 
by Eq. ([1]), we write the magnetic energy per lattice site 
for an arbitrary 2D propagation vector k = (k x , k y ): 

E = - (J cos k x + J\ cos k y + J 2 cos(2fc y )) , (6) 
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FIG. 5. (Color online) Left panel: the pitch angle (7) and 
the sublattice magnetization (m) calculated by the CCM for 
J2 = -§|Ji|. The shaded bar shows the coupling regime of 
Cu2Ge04. Note that the quantum pitch angles 7 qu for the 
LSUBn approximations with n = 4, 6, and 8 almost coincide. 
Therefore, the shown LSUB8 curve represents effectively the 
limit n — > 00. Right panel: spin-spin correlation (SoSr), 
R = (x, y) along the Jj — J2 chains (fe'-axis) for a finite system 
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where the unit cell of the spin lattice is used<^ The en- 
ergy minimum is found at k = (tt, arccos(— jj-)) that 
corresponds to the AFM order along a' and spiral order 
along b' . The classical pitch angle is 7 = arccos(— 4j-) = 
79.19° and does not depend on J, hence the a' and b' 
directions of the spin lattice are fully decoupled. The or- 
dering along a' is controlled by J, whereas the ordering 
along b' is controlled by the competing couplings J\ and 
J 2- To check the validity of this result for the quantum 
case, we use the CCM method and Lanczos diagonaliza- 
tion. 

In CCM, we reduce our exchange couplings (Tabled to 
Ji = -1, -h = |, and vary J (in Cu 2 Ge0 4 , J = f ). The 
CCM results for 7 and to as a function of J are shown 
in Fig. [S] In contrast to the classical pitch angle 7 c i, the 
pitch angle of the quantum system (7 qu ) slightly depends 
on J. In Cu2Ge04, we find 7 qu = 83.9°, which is about 
6% larger than the classical angle, but about 5% smaller 
than the quantum pitch angle for the isolated chain, i.e., 
at J = 0. The coupling J affects the dimensionality of 
the system and has a stronger effect on the sublattice 
magnetization (see Fig. [S]). The extrapolated value Toco 
has a maximum at J ~ — 1.17Ji. The calculated ex- 
change couplings in Cu2Ge04 lead to TOqo ~ 0.310 (i.e., 
0.62 n B ). 

The CCM results are confirmed by the Lanczos diag- 
onalization data for the spin-spin correlation functions 
(SoSr) shown in the right panel of Fig. [5J We use a fi- 
nite lattice comprising four 8-spin J\ — J2 chains coupled 
by J. The correlations between nearest neighbors within 
the frustrated J1 — J2 chains [R = (0, 1)] are close to zero, 
whereas the second-neighbor correlations [R = (0, 2)] are 
AFM. This conforms to the spiral ordering with the pitch 
angle close to 90° (neighboring spins are nearly orthogo- 
nal). The correlations between the structural chains [at 
R = (1, y)] follow the intra-chain correlations, yet show- 
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ing the opposite sign. Thus, the ordering along a' is 
AFM. 



D. Long-range order 

The 2D model given by Eq. ((I) is ordered at zero tem- 
perature only. To account for the actual long-range order 
in Cu 2 Ge04 below T/v = 33 K, the interlayer coupling 
J c should be considered. The FM coupling J c is compat- 
ible with Ji, yet competing with J and J2. Assuming a 
similar ground state with the 2D propagation vector, we 
find that J c modifies the energy in Eq. ([6]) by 

A£ 3D = ^cos^cos^. (7) 

Using DFT estimates of individual exchange couplings 
(Table HI, we arrive at the classical pitch angle modified 
by 0.15 %: 79.07° vs. 79.19° for the purely 2D model. 
The classical energy per lattice site is reduced by 0.5 % 
(about 0.6 K). This simplified analysis shows that the 
interlayer coupling J c is capable of stabilizing the 3D 
order in CU2GCO4. However, the classical model does 
not reflect all the features of the real quantum model 
that, unfortunately, remains unfeasible for an accurate 
numerical study. In particular, the ordering temperature 
Tjv can not be determined with sufficient accuracy. 

IV. DISCUSSION AND SUMMARY 

Although not obvious at first glance, the microscopic 
magnetic model of Cu2Ge04 can be deduced from sim- 
ple qualitative arguments. While the naive geometrical 
analysis of the crystal structure suggests a 3D pyrochlorc- 
lattice magnetism, a closer look at the crystal structure 
identifies 2D features. In most of the Cu +2 oxides, elec- 
tronic structure and magnetism are controlled by the ar- 
rangement of CUO4 plaquettes, which are the basic struc- 
tural entities. Chains of edge-sharding plaquettes give 
rise to frustrated 3\ — J 2 spin chains,— ~— yet the coplanar 
arrangement of the plaquettes in the neighboring struc- 
tural chains induces a strong AFM coupling along a'. 
Overall, we find magnetic layers in the ab plane, along 
with a weak FM interlayer coupling J c . By combining 
the frustration along b' with the strong unfrustrated ex- 
change along a', Cu2Ge04 expands the family of cupratcs 
featuring frustrated spin chains. 

The dearth of the experimental data and the complex- 
ity of the 2D frustrated J — J\ — J 2 lattice restrict the 
opportunities for an experimental verification of our mi- 
croscopic model. Nevertheless, the tangible success of 
DFT in unraveling complex spin lattices for a range of 
transition-metal compound a 27 i 35 ' 51 i 52 is a solid justifica- 
tion of our results. The 2D nature of the system and the 
frustrated couplings along b' are confirmed by qualita- 
tive arguments and by a reference to similar Cu +2 com- 
pounds (see Sec. |n| . Further on, numerical estimates of 



individual exchanges conform to the experimental mag- 
netic susceptibility (Sec. IIII A|) . An ultimate test of the 
proposed model requires a study of the ground state by 
neutron or resonant x-ray scattering. Presently, we note 
that our model does predict the long-range magnetic or- 
der, in contrast to the strongly anisotropic pyrochlore 
lattice that might have a quantum-disordered valence- 
bond-solid2^ or a gapless spin-liquid^ ground state. 

The experimental data for Cu 2 Ge04 and accurate the- 
oretical results for the ground state disclose the basic 
features of our model. At high temperatures, thermody- 
namic properties are guided by the uniform spin chains 
along a' . The apparent spin-chain behavior is likely re- 
lated to the partial cancellation of FM J x and AFM J 2 in 
the second-order term for the susceptibility (see Eq. Q). 
A further evidence is the perfect fit of the experimental 
magnetic susceptibility down to 70 K (T/J ~ 0.5). The 
frustrated couplings along b' come into play at lower tem- 
peratures, and essentially determine the ground state. 
The leading exchange J drives AFM ordering along a'. 
The ordering along b' has to satisfy the frustrated cou- 
plings Ji and J 2 and is, therefore, a spiral, similar to a 
single Ji — J 2 frustrated spin chain. The interlayer cou- 
pling J c should stabilize the long-range order up to Xjv 
without changing the basic features of the ground state: 
the collinear AFM order along a' and the spiral order 
along b'. 

The highly accurate CCM approach provides reliable 
information on the ground state of the 2D system. We 
find the pitch angle of 7 ~ 84° and the sublattice magne- 
tization close to 0.62 Both 7 and m are renormalized 
with respect to the classical values, and suggest strong 
quantum fluctuations in the system. Enhanced quantum 
fluctuations should be ascribed to the reduced dimen- 
sionality and frustration. The magnetic ordering temper- 
ature T/v / J — 0.25 is also suggestive of strong quantum 
fluctuations. For example, a quasi-2D system of square 
lattices with a weak interlayer coupling J±/J = 0.01 
(compare to \J C \/J — 0.01) orders at a higher temper- 
ature of Tn / J ~ 0.33 (Ref. [54|). The quantum effects in 
the system could be further probed by an experimental 
study of the ground state. 

An interesting feature of the spiral magnetic order is 
the possible emergence of electric polarization strongly 
coupled to the magnetism^ The direction of the electric 
polarization depends on the twisting direction of the spi- 
ral. In Cu2Ge04, the AFM coupling J leads to opposite 
twisting directions in the neighboring spirals, hence the 
polarization is canceled. However, the proposed antifer- 
roelectricity of Cu 2 Ge04 does not preclude the strong 
magnetoelectric coupling, and should stimulate further 
experimental investigation of the compound. We also 
note that the family of 2D frustrated materials repre- 
senting the J — J\ — J2 model can be further expanded 
by CuNCN lying in the limit of J » \Ji\, J 2 M 

In summary, we have shown that the electronic struc- 
ture of CU2GCO4 contradicts the previous, empirical- 
based spin model of the anisotropic pyrochlore lattice. 



7 



The comprehensive computational study discloses the 
quasi-2D nature of this compound and suggests an orig- 
inal 2D spin model comprising frustrated and uniform 
spin chains along the two dimensions. Theoretical results 
for this model show the robust nature of the spiral ground 
state that is subject to strong quantum effects, evidenced 
by the reduced sublattice magnetization of 0.62 [Ib and 
the renormalized pitch angle of about 84°. 
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